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Abstract 

We study direct searches of additional Higgs bosons in multi-top-quarks events at the LHC Run- 
H, its luminosity upgraded version with 3000 fb“^, and the International Linear Collider (ILC) with 
the collision energy of 1 TeV. Additional Higgs bosons are predicted in all kinds of extended Higgs 
sectors, and their detection at collider experiments is a clear signature of the physics beyond the 
standard model. We consider two Higgs doublet models with the discrete symmetry as benchmark 
models. If these additional Higgs bosons are heavy enough, the decay modes including top quarks 
can be dominant, and the searches in multi-top-quarks events become an important probe of the 
Higgs sector. We evaluate the discovery reach in the parameter space of the model, and find that 
there are parameter regions where the searches at the LHC with 3000 fb“^ cannot survey, but the 
searches at the ILC 1 TeV run can. The combination of direct searches at the LHC and the ILC 
is useful to explore extended Higgs sectors. 
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I. INTRODUCTION 


A Higgs boson (h) was discovered at the LHC Run-I in 2012 After the discovery, fur¬ 

ther precision measurements have revealed its properties, which seem to be quite consistent 
with those of the Higgs boson in the standard model (SM) within current experimental er¬ 
rors By measuring the mass of the Higgs boson very precisely, rrih = 125.09±0.24 GeV 
at the LHC 0], the parameters in the Higgs potential in the SM have been determined ex¬ 
perimentally. However, the whole structure of the Higgs sector is still unclear, since one 
can consider models with extra scalar helds, like the Higgs sector in the supersymmetric 
models, satisfying all the available constraints including the properties of the discovered h. 
Moreover, in fact, numbers of models with extra scalar helds have been proposed to solve the 
problems which cannot be explained within the SM, such as neutrino masses, dark matter, 
baryon asymmetry of the universe, cosmic inhation, etc. Determining the structure of the 
Higgs sector experimentally becomes an important foothold in constructing the new theory 
above the electroweak scale. 


For this purpose, further precision measurements are required to pin down small devia¬ 
tions from the SM in the properties of h P, This shall be performed at the LHC Run-H 
with increased energy and more accumulated luminosity, and also at the future lepton col¬ 
liders, such as the International Linear Collider (ILC) 0,0, the Future Circular Collider 
(FCC) 10, and the Compact Linear Collider (CLIC) [0. By measuring deviations from 
the SM in the properties of h, one can probe an extended Higgs sector with rather large 
masses of additional Higgs bosons so that effects by the extended part in the Higgs sector 


tend to decouple from the properties of the SM-like h [14, 15 


Another approach to elucidate the whole structure of the Higgs sector is to directly search 
for additional Higgs bosons, since their existence must be a clear evidence of an extended 
Higgs sector. Direct searches of additional Higgs bosons have been performed at the LHC 


Run-I in various decay modes |0, and the limits on the mass of additional Higgs bosons 
and their coupling strength have been investigated. On the prospect of the direct searches at 
future colliders, the LHC Run-H has an advantage since it is the energy frontier experiment 
which is good to produce heavier particles. However, on the other hand, direct searches at 
the ILC have a different advantage, that is, a parameter regions with relatively small cross 
section can be probed due to the clean environment of lepton collider experiments, although 
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the mass reach for additional Higgs bosons is relatively limited. Thus, the searches at the 
LHC Run-H or its luminosity upgraded version with 3000 fb“^ (LHC 3000 fb“^) and the ILC 
can be complementary to servey the wide parameter regions in extended Higgs sectors [I^ . 


In this paper, we study multi-top-quarks events as distinct signals of production of ad¬ 
ditional Higgs bosons. As benchmark models of extended Higgs sectors, we consider two 
Higgs doublet models (2HDMs).Q In the 2HDM, there are CP-even H, CP-odd A, and a 
pair of charged in addition to the SM-like h. For neutral H and A, the decays into a 
top-quark pair open if their masses are larger than about 350 GeV. Because of the large mass 
of the top quark, it is quite possible that this decay mode dominates the branching fraction 
at least in certain regions in the parameter space. Thus, the multi-top-quarks events can 
be an attractive signal of the heavy Higgs bosons. Such multi-top-quarks events have been 
studied at the LHC as a signal of new particle in various models [^-^|. A large production 
rate of multi-top-quarks events is predicted at the LHC for the models with colored new 
particles such as gluinos, color-octet scalars, etc., but that for non-colored particles such as 
heavy additional Higgs bosons is limited. At the ILC, four top-quarks production can be 
considered as a signal of H and A, through —?• HA and f fH/A [|n 


36 


To produce a pair of H and A which decay into top-quark pairs, both the masses of H and 
A are required to be larger than about 350 GeV, and the collision energy should be higher 
than about 700 GeV. Such an experiment can be realized at the ILG with ^/s = 1 TeV. 
Up to our knowledge, there has been no dedicated study on this process at lepton colliders. 
Therefore, in this paper, we aim to present a detailed analysis on this process including the 
hadron-level simulation with jet clustering, flavor tagging, detector acceptance and momen¬ 
tum resolution effects. We hnd that the four top-quarks events can be detected by simple 
kinematical cuts, and thus be useful to survey the parameter regions in the 2HDM. We note 
that the signal of heavy charged Higgs boson can be —)■ tb{tb), and its observability has 


been studied in gb —?■ process at hadron colliders |^, and also at lepton colliders 

in e+e- ^ H+H- and e+e" ^ H^ff 

The paper is organized as follows. In Sec. |1|, we briefly introduce the 2HDM with Z 2 
symmetry considering the four types of Yukawa interactions. In Sec. B we present an 
analysis for the search prospect of additional Higgs bosons in multi-top-quarks events at 
the LHG. In Sec. |V|, we study the four top-quarks events at the ILG by performing the 
^ For the review of the 2HDMs, see, e.g., Refs. [Q 
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TABLE I. Four possible Z 2 charge assignments to the scalar and fermion fields. 


Monte-Carlo simulation for the signal and background processes at the hadron level with 
detector effects. By using the simulation analysis, we evaluate the discovery potential of the 
neutral Higgs bosons at the ILC in the four top-quarks events in the parameter space in the 
2HDM with four types of Yukawa interactions. The obtained discovery reaches at the LHC 
and at the ILC are compared. Sec. 0 is devoted to discussions for further investigation and 


future prospects. Finally, we draw a conclusion in Sec. ^ 


II. TWO HIGGS DOUBLET MODEL 


In this section, we briefly introduce the model we consider, namely the 2HDM with the 
softly-broken symmetry. We introduce two isospin doublet scalar helds, $1 and <1)2, which 
transform as $1 —)■ -|-<I)i and ^2 —<1)2 under the Z 2 transformation. For the SM fermions, 

there are four kinds of Z 2 parity assignment [42-Q, as listed in Table |. We denote the four 
types of Yukawa interactions as Type-I, Type-H, Type-X and Type-Y The Yukawa 

interaction to the SM fermions in each flavor is allowed for only one Higgs held, <Fi or <^ 2 , 
to make each interaction term Z 2 invariant. It forbids the havor changing neutral currents 


at the tree-level [45 , which are severely constrained by experimental observations. 


After solving the condition for the electroweak symmetry breaking and diagonalizing the 
mass matrices in the Higgs sector assuming CP-invariance, there are hve physics scalar 
states; namely, two CP-even neutral Higgs boson h and P, one CP-odd neutral Higgs 
boson A, and a pair of charged Higgs bosons H^. The lighter CP-even h can be identihed 
as the SM-like Higgs boson which was observed at the LHC. Yukawa interactions are given 
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TABLE II. The scaling factors 
sm0 for 9 = a, (3. 


for the four types of Yukawa interactions | ^ ]. ce 


cos 9, and sq = 


in terms of these physical scalars as 


-£v«k.w.= Y, + 

I V V V 


f=u,d,i 


[tti^CaPl + dH+ + + H.c. \ , (1) 


where the scaling factor with cj) = h, H, A and / = u,d,i can be found in Table |l[ The 
scaling factor is a function of a and [3, the mixing angles in the neutral CP-even component 
and CP-odd component, respectively. The mixing angle (3 also satishes tan f3 = ^ 2 /'^!, where 
Vi and V 2 are the vacuum expectation values of the two doublet helds. 

The gauge coupling of h is given by glvv^ — 9h^v sin(/5 — a) and that of H is given by 
— 5'fcvv cos(/5 — a). Theoretically, a deviation of sin(/3 — a) from unity is constrained 
by the arguments of perturbative unitarity and vacuum stability |^9|-p|. If a soft- 

breaking scale of the Z2 symmetry M is larger than the electroweak scale, M S> n, only 


small value of 1 — sin(/9 — a) is allowed by these constraints . The limit of sin(/5 — a) —)■ 1 
is called the SM-like limit, where h has the same coupling constants to the gauge bosons and 
also to the SM fermions as the SM Higgs boson. On the other hand, Yukawa interactions of 
P, A and to the SM fermions do not vanish in this limit, and the coupling strength for 
each vertex depends on the type of Yukawa interactions and tan/?. Thus, the variety of the 
type of Yukawa interactions with different tan (3 dependences leads to rich phenomenology 
for the additional Higgs bosons. 

We focus on the interactions of additional Higgs bosons to top quarks. For any type of 
Yukawa interactions, the Yukawa coupling constants to top quarks are enhanced by cot (3 for 
small tan (3 regions. Therefore, for larger masses of additional Higgs bosons where the decay 
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modes into top quarks are open, the branching ratio to top quarks becomes dominant for 
small tan (3 regions. The hgures for the branching ratio of additional Higgs bosons can be 
found, e.g., in Refs. BEl- For larger tan /3, the dominant branching ratio is replaced by 
the other fermionic mode, bb for Type-II and Type-Y, for Type-X Yukawa interactions. 
For Type-I, since the tan /3 dependence is common for all fermions, the dominance of ti decay 
mode is true for any value of tau/d. 


III. MULTI-TOP-QUARKS PRODUCTION AT THE LHC 


In this section, we study the four top-quarks production through the production of addi¬ 
tional Higgs boson(s) at the LHC. The largest contribution comes from the top-quark pair 
associated production process, 

pp —>■ tiH{tiA) —)■ titi, (2) 

since it emerges via the strong interaction. At the tree level, there are two subprocesses 
in this process. One is gg —)■ tiH{ttA) and the other is qq —)■ ttH{ttA). Because H{A) is 
radiated off from the top quarks, the cross section is proportional to the square of {yf) 
which is proportional to cot (3 in the SM-like limit. Therefore, the cross section is large for 
smaller tan/?. 

The four top-quarks production through the pair production of H and A, 

pp —>■ HA —>■ tttt, (3) 


is described by the quark anti-quark annihilation process, qq ^ Z* ^ HA aX the tree level. 
Since the HA production cross section does not depend on tan/?, the cross section of the 
hnal four top-quarks production depends on tan (3 through the branching ratios of H and A 
into the top-quark pair, 

There are also three top-quarks production processes via the associated production of 
H^ and H{A), which subsequently decay into tb(tb) and ti, respectively, 

pp —)■ tbti/ibti. (4) 


At the tree level, this production process is described by the W boson mediated diagram |53 
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qq' ^ IF* ^ H^H{H^A). 


( 5 ) 
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We estimate the cross sections for these processes at the LHC at the tree level. As an 
example, we take the 2HDM with Type-II Ynkawa interactions. For simplicity, we take a 
common mass for all the additional Higgs bosons. Onr calcnlation is performed with the 
nse of analytic eqnations in Ref. [Q and the nnmerical codes generated by MadGraphS 


To calcnlate the branching fractions oi H ^ ti and A —)• ti, the off-shell effect of top 
qnarks is inclnded. To estimate the hadronic cross section, the CTEQ6L parton distribntion 
fnnctions (PDFs) are nsed with setting the scale of PDFs to /i^ = for HA 

and H^H -|- H^A prodnction processes, and to /ip = m* -|- m$/2 for the ttH (A) prodnction 
process In Fig. |I], the resnlts are plotted as a fnnction of the mass for the LHC 8 TeV 
(left) and the LHC 14 TeV (right). The cross sections for ttH -|- ttA, HA and H^H + H^A 
prodnctions are plotted in bine, black and red dotted lines, respectively. For the ttH + tiA 
prodnction process, the resnlts for tan/3 = 1 and 3 are plotted. In addition, the resnlting 
cross sections for the mnlti-top-quark prodnction processes are plotted in solid (dashed) lines 
with the same color, assnming that tan/3 = 1 and 3. The cross sections at 13 TeV collision 
energy are abont 70-80% of those at 14 TeV. 

For each process, the largest cross section is realized for the mass of additional Higgs 
bosons at aronnd 350 GeV. Below this valne, the branching ratio into ti is snppressed 
becanse one of the top qnarks is forced to be off-shell, while above that valne the prodnction 
cross sections of additional Higgs bosons get decreased. For tan/3 = 1, the cross section of 
the fonr top-qnarks prodnction can be at most 6 fb (50 fb) for the LHC 8 TeV (14 TeV). 
We note that there has been already an experimental npper limit for the cross section of 
fonr top-qnarks prodnction at the LHC 8 TeV |^, The CMS Collaboration has set 
the limit to < 32 fb at the 95% CL (conhdence level) |^, by observing the lepton pins 
mnlti-jets events. However, the limit is not tight enongh to constrain the parameter regions 
in the 2HDM. We also note that the SM prediction to the fonr top-qnarks prodnction at 


the LHC 8 TeV is abont 1 fb For the three top-qnarks prodnction, since the expected 
cross section is smaller than that of fonr top-qnarks prodnction while the signatnres mix np 
with that of fonr top-qnarks prodnction, the detection may be more challenging. 

We now stndy the prospect of measnring fonr top-qnarks prodnction as a signal of the 
prodnction of additional Higgs bosons at the fntnre LHC rnn. For the fonr top-qnarks 
prodnction within the SM contribntion, the total cross section and its nncertainty are esti¬ 
mated to be (Tsm = 15 fb and (Jctsm = 4 fb, respectively |Q. To deal with the backgronnd 
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FIG. 1. Cross sections for multi-top-quarks production at the LHC 8 TeV (left) and 14 TeV (right) 
for the Type-II 2HDM. four top-quarks production from ttH + ttA production, HA production 
and three top-quarks production from H-j- AH^ production are shown with tan /3 = 1 (solid 
lines) and 3 (dashed lines). 


processes, we follow the analysis in Ref. where the selection cnts to extract the fonr 
top-qnarks events ont of the background events are demonstrated by simulation analysis. 
In their analysis, the background rate of B = 7.2 fb after selection cuts is obtained with 
the signal efficiency of e = 0.03. By taking into account the statistical and systematical 
uncertainties for the signal, SM and background processes, the accuracy of measuring the 
signal cross section as can be estimated as 


S(Ts I (crs + ctsm) e + 5 5aly^e^ -k {dBf 

Ty-y—;pz— + —;p—■ 

where 5B denotes the systematic uncertainty of the background rate. We take 5B = O.OSR, 
which may be achieved at the later stage of the LHC experiment. By solving Eq. (0) , we 
obtain that as has to be larger than 25 fb (63 fb) to achieve Sas/cTs < 0.5 (0.2) with the 
integrated luminosity of C = 300 fb“^. In our setup, the total uncertainty is dominated 
by the systematic uncertainty from the background. To reduce the statistical uncertainty 
smaller than the systematical one, one needs only more than 10 fb“^ of the data. Thus, the 
accuracy will not be improved by accumulating the integrated luminosity up to 3000 fb“^, 
but is limited by the systematical errors. Increased integrated luminosity is, on the other 
hand, useful to reduce the systematical uncertainty in the backgrounds. If we change our 
input by R = 7.2 fb ^ 3.6 fb or e = 0.03 —)■ 0.06, the resulting values are modihed to 
14 fb (36 fb) for Sas/crs < 0.5 (0.2), respectively. In the ideal case of 5B = 0, the parameter 
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FIG. 2. Contour plot for the discovery reach of four top-quarks events at the LHC Run-II in the 
2HDM. Upper limits on tan (3 for the discovery regions at the 2cj [5 cj] CL are shown in dashed 
[solid] lines. 


regions with cr^ > 8 fb (20 fb) can be observed at the 2cr (5cr) CL. To obtain better sensitivity, 
a considerable amonnt of efforts to rednce the systematical nncertainties is reqnired. 

In Fig. 1^, we show the parameter regions in the 2HDM with fonr types of Ynkawa inter¬ 
actions where the additional Higgs bosons contribntion can be detected at the LHC in the 
fonr top-qnarks events at the 2a and 5a CL for the above setnp. It tnrned ont that the 
dependence on the type of Ynkawa interactions is small, and only tan/3 < 1.5 can be probed 
at the 2a CL at most. Since these regions are constrained by flavor experiments [^, |6^ , 
the LHC searches in 4 top-qnarks events may not have signihcant impact on exploring the 
parameter regions in the 2HDMs. 
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IV. MULTI-TOP-QUARKS PRODUCTION AT THE ILC 


In this section, we consider the four top-quarks production at the ILC. In the 2IIDM, the 
four top-quarks hnal state is generated via the pair and single production of H and/or A. 
For s/s > mriH + pair production of H and A, 

e+e" ^ HA, (7) 

is kinematically possible and its cross section can be sizable. The HA pair production cross 
section does not depend on tan/? at the tree level. Thus, the four top-quarks production 
rate depends on tan /? only through the decay branching ratio of H and A. 

On the other hand, for ^/s < mn + rnA, the pair production is kinematically forbidden, 
but the single production process 


e’''e -^ttH(ttA), 


( 8 ) 


can still open as long as ^/s > 2mt+mH{A)- The cross section of this process can be increased 
by the enhanced Yukawa coupling of H and A to the top quarks. Followed by the decays of 
H and A into ti, four top-quarks events occur from these processes. 

Through the decay of top quarks, the signature of four top-quarks production can be ob¬ 
served as all-hadronic, single lepton plus jets plus missing momentum, dilepton plus jets plus 
missing, trilepton plus jets plus missing, tetralepton plus jets plus missing channels. Among 
dilepton plus jets plus missing channels, there are same-sign and opposite-sign dilepton hnal 
states, where the former is expected to have small backgrounds. The branching fractions 
for these channels are listed in Table |T|. In the SM, the leading production mechanism of 
four top-quarks events is —)■ ttg* —)■ titi via QCD interactions, thus the cross section 

is The next-to-leading production mechanism is full electroweak process, 0{a^). 

When we consider the hnal states including the decay of top quarks and their detection at 
real experiments, there enter reducible backgrounds via —)■ ti, tibb, etc. We 

estimate the contributions from these processes in the following simulation analysis. 

First, we describe our simulation analysis for the detection of the 4 top-quarks events 
at the ILC with ^/s = 1 TeV. We use MadGraphS [Q and PythiaG for generation 
of signal and background events, with Tauola m for tan lepton decays. We evaluate 
the signal process by including both the single and pair production amplitudes coherently. 
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Decay modes 

Final states 

TZ (with r’s) 

TZ (without r’s) 

all-hadron 

+ 8j 

(ir-0-2 

(ir-0-2 

single lepton -|- jets 

+ Ajb + 6j + 1^ 

(|)^W4c.0.4 

(i)^ •1-4-0.26 

S.S. dilepton -|- jets 

+ Ajb + 4:j -|- UE 

(i)'-0'-2-o.i 

(I)'-(i)'-2-0.04 

O.S. dilepton -|- jets 

+ Ajb + 4j -|- EE 


(i)'-(i)'-4-0.09 

trilepton -|- jets 

£^£^£^ + Ajb -f 2j + EEE 

l-Qf-^-o.i 

|. (1)^4-0.03 

tetralepton -|- jets 

£'^£'^£~£~ -|- Ajb + EEEE 

(ir-o-oi 

(1)"^ ~ 2.4 X 10"3 


TABLE III. Branching fractions of the four top-quarks signature in the decays of top quarks. 
Branching fractions where includes or not are listed. 

The SM contribution to the e’''e“ —?• tttt process is estimated to be very small, giving 
o'tot = 3.8 X 10“^ fb at y/s = 1 TeV. Thus, for the integrated luminosity of £ = 1 ab“^, only 
a few events are expected to be produced in the SM. The background processes considered 
in our analysis are 


e’'"e“ —)■ ti, 

( 9 ) 

e'^e~ —)■ ttbb, 

( 10 ) 

e'^e~ tti^i~. 

( 11 ) 

e+e" ^ ttW+W-. 

( 12 ) 


The second process includes e+e“ bb), e+e“ —)■ bb), e~^e~ —)■ ttZ/ 7 *(—)■ bb), 

e~^e~ —)■ tbW*{^ tb), and e+e“ —>■ W~^*W~* ttbb. in the third process mean the sum 
of e^, /i^ and r^. The other background processes are negligible. 

To analyze the gererated events, we follow the designed performance of the ILC detec¬ 
tors E|. We take all detectable particles whose pseudo rapidity satishes l^l < 1.5. For 
charged particles, we further require their transverse momentum satishes pt > 0.3 GeV. 
Four momenta of those particles are smeared by using the Gaussian distribution with 
^pt/Pt = + 10“^ for charged particles, (Te/E = 0.4/\/^ -|- 0.02 for neutral hadrons, 

and cfe!E = 0.15 /\fE -|- 0.01 for photons, where pt and E are given in an unit of GeV. 

Among those particles, we select isolated leptons, e and p, whose energy satishes i^cone < 
^^/Q{Ee — 15) where Econe and Ei are given in an unit of GeV. Econe is the summation of 
energies of particles inside the cone around the lepton dehned as cos^cone > 0.98 except the 
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FIG. 3. Kinematical distributions of signal and background events in Thrust (left), (middle), 
and Nt = 2iViep + A(jet (right) variables with loose 6-tagging criteria. Signal events are displayed 
for mH,A = 400 GeV with a normalization of crtot = 2.7 fb which corresponds to tan/3 = 1 in 
Type-II 2HDM. 


lepton itself 

After removing the isolated leptons from the list of particles, we perform a jet clustering 
by using the Durham algorithm with hxed Dcut = 5 x 10“^ with the help of Fast jet ^ . 


Thus, the number of jets in an event is flexible according to the event structure. For each 
clustered jet, we perform a flavour tagging. If a jet contains only photons, it is tagged as 
a photon-jet. i3-tagging is performed stochastically by using the decay history of particles 
which is available in the Monte-Carlo simulation. If a jet contains S-hadrons (D-hadrons) 
in the decay history of constituent particles, we tag it as a b-jet randomly with a probability 
of 80% (10%). A probability of mis-tagging a jet which does not contain B or D-hadrons 
as a b-jet is set to be 3%. These probabilities correspond to loose tagging criteria given in 
Ref. [^. We found that this loose criteria works better than rather tight criteria to collect 
more signal events in the circumstances of small backgrounds. In addition, a jet is tagged 
as a tau-jet, if it contains 1 or 3 charged tracks and satishes Econe/Ejet > 0.95 where Econe 
is the summation of energies inside the small cone around a direction of jet three momenta 
with R = 0.15. The other jets are assumed to be light-jets. The number of leptons in an 
event is counted as Wep = ° + Wj, where is the number of isolated e (/i) and 

Nt-. is the number of tau-jets. The number of jets in an event is counted as Ajet = Nij N^j, 
where Nij^pj) is the number of light-jets (b-jets). 
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Processes 

Cross sections 

Accumulated efficiencies 


o-tot [fb] 

T < 0.77 

Vfc. > 3 

2Ne + Nj > 10 

e''“e (— )• HA) ^ tttt [Typ 

3-II, sin(/3 - 

- a) = 1, tan /3 = 1] 

mH,A = 360 GeV 

4.3 

71% 

58% 

34% 

400 GeV 

2.7 

92% 

74% 

43% 

440 GeV 

1.3 

96% 

79% 

47% 

480 GeV 

0.30 

96% 

77% 

46% 

500 GeV 

7.5 X 10-2 

95% 

77% 

45% 

520 GeV 

3.2 X 10-2 

95% 

77% 

45% 

560 GeV 

1.0 X 10-2 

93% 

76% 

44% 

e+e” -A- tt 

166. 

3.6 X 10-3 

5.6 X 10-“^ 

2.0 X 10-6 

e~^e~ -A ttbb 

5.0 

19% 

14% 

0.66% 

e+e“ ^ ttNi~ 

0.76 

23% 

3.8% 

1.4% 

e+e- -A tiW+W- 

0.14 

55% 

15% 

3.0% 

e'''e“ —)■ tttt (SM) 

3.8 X 10-3 

93% 

74% 

41% 


TABLE IV. The total cross section and the accumulated efficiencies by kinematical cuts for the 
signal and background processes at the ILC with ^/s = 1 TeV. Signal cross sections and efficiencies 
are calculated for mH,A = 360 GeV to 560 GeV with tan /3 = 1 for Type-II as a reference. 


To extract the signal events ont of the SM backgronnd, we impose following selection cnts; 
(1) small thrnst, T < 0.77, (2) Ni,j > 3, and (3) large '^hard multiplicit'tj\ Nt = 2iViep+Ajet > 
10. In Fig. D, we show distribntions of signal and backgronnd events in thrnst, N^j and Nt 
in the left, middle and right panels, respectively. Backgronnds from ti, and ttbb 

productions are shown in blue, cyan, and red histograms, respectively. Contributions from 
the tiW^W~ process and the SM four top-quarks production are omitted since these are 
small enough. For the reference, signal events with mH,A = 400 GeV are shown on top of 
them where the normalization is adjusted to the case of tan/3 = 1 in the Type-II 2HDM. 
We hnd that large tt backgrounds are excluded by the thrust cut, and a large amount of 
tt, contributions is excluded by the Nbj cut. Furthermore, all the backgrounds are 

suppressed by applying the cut on Nt- 


The background reduction and signal detection efficiencies are summarized in Table IV 
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where the background process Eq. (^) and the SM four top-quarks production processes 
are also included for the reference. In our simulation, the background reduction rates by the 
above three cuts are 0(10“®) for tt, 0.66% for ttbb, 1.4% for and 3% for ttW~^W~. 

With the integrated luminosity of £ = 1 ab“^, only around 47.8 events are expected to be 
observed. From the SM four top-quarks production, we expect 1.6 events. On the other 
hand, for the signal process of additional Higgs bosons, around 34% to 47% of events are 
remained after these cuts, depending on the mass. 

By using the signal detection efficiency es with the expected number of background rate 
B = 49.4 ab at £ = 1 ab“^, we can estimate the minimum value of the total cross section 
for the signal process to be identified in a certain accuracy. We take into account only 
the statistical uncertainties due to the presence of backgrounds, since at lepton colliders 
systematical uncertainties can be expected to be well under control. Thus, the uncertainty 
of observing the total cross section of signal events is estimated as 

— = ( 13 ) 

We find that the signal process can be detected at the 2a (ha) CL if the total cross section 
is above 0.034-0.048 fb (0.11-0.15 fb), depending on the mass of additional Higgs bosons 
through the signal detection efficiency. 

In Fig. we evaluate the parameter regions in the tan/?) plane in the 2HDMs 

where the four top-quarks events can be detected at the ILC. Solid (dashed) contours show 
the regions where the signal can be detected at the 2a (5a) CL. For Type-H and Type-Y 
2HDM, tan/? up to around 7 (6) can be probed for the mass of additional Higgs bosons 
up to mH,A ^ 500 GeV. Above mH,A — 500 GeV, the signal can be produced by using the 
single Higgs production processes, Eq. (H), but only visible for small tan/? below unity. For 
Type-X 2HDM, tan/? up to around 15 (12) can be probed for mu ,a ^ 500 GeV, because of 
the fact that the H/A -A- ti decay branching rate remains to be dominant up to larger tan /? 
than that in Type-H and Y. For Type-I 2HDM, since the H/A ^ ti decay mode dominates 
for any value of tan /?, the four top-quarks events can be observed for any value of tan /?, as 
long as mH,A ^ 500 GeV. 

We note that, however, the result for Type-I is due to the fact that we take the SM- 
like limit, sin(/? — a) = 1, where no 77 —>■ WW, ZZ, hh or A —>■ Zh decay is induced. 
For example, if we take sin(/? — a) = 0.99, these decay modes become non-zero, and even 
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FIG. 4. Contour plot of the discovery reach of the four top-quarks production in the 2HDMs at 
the ILC ^/s = 1 TeV with C = 1 ab“^ data. Upper limits on tan /3 for the discovery regions at the 
2a [5cj] CL are shown in solid [dashed] lines. The similar results but for sin(/3 — a) = 0.99 are also 
shown in blue lines. 


dominant for larger tan/3 |T^. The value of sin(/3 — a) = 0.99 can be detected by the 
precision measurement of hVV coupling constants at the ILC with ^/s = 500 GeV and 
C = 1.6 ab“^ at the 5a CL [§]. In this case, the four top-quarks events can be observed 
up to about tan /3 ~ 8 (5) at the 2a (5(j) CL for the four types of Yukawa models. The 
discovery reaches for sin(/3 — a) = 0.99 are also shown in each hgure in Fig. ^ in blue lines. 


We now compare the impact of the searches for multi-top-quarks events at the LHC 
and the ILC on exploring the parameter regions in the 2HDM. At the LHC Run-H, by the 
searches for four top-quarks events, the parameter regions with tan/3 < 1.5 can be surveyed 
in all types of Yukawa interactions. Because the experimental uncertainties are dominated 
by the systematical ones, the results will not be improved at the LHC 3000 fb“^. On the 
other hand, at the ILC 1 TeV, we hnd that the parameter regions with larger tan /3 can be 
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surveyed as long as mH,A ^ 500 GeV. For Type-II and Type-Y (Type-X), the parameter 
regions with tan/? ~ 7 (15) can be explored. For Type-I, because the decay mode into 
a top-quark pair is dominant for any value of tan /3, detection of four top-quarks events 
is anticipated in any value of tan/?. This result is modihed to tan/? < 9 in the case of 
sin(/? — a) = 0.99. Therefore, the searches for four top-quarks events can explore wider 
parameter regions at the ILC 1 TeV than at the LHC 3000 fb“^. 


V. DISCUSSIONS 


In this section, we discuss further investigation on the multi-top-quarks events at the 
LHC and the ILC, regarding the mass determination of the additional Higgs bosons and the 
discrimination of the type of Yukawa interactions. In addition, we give some comments on 
the searches for multi-top-quarks events at future multi-TeV lepton collider, CLIC. 

First, we discuss the mass reconstruction of additional Higgs bosons in the decay mode 
into a top-quark pair. At the LHC, even if some excess in the events for the lepton plus 
multi-jets channel is observed, the mass reconstruction of additional Higgs boson is difficult, 
since one has to form an invariant mass of six jets in the pretense of more jets which come 
from the top quarks and the ISR/FSR jets. The combinatorial uncertainty and the jet energy 
resolution prevent us from constructing a clear peak in the invariant mass distribution. At 
the ILC, the mass reconstruction by kinematical methods is still difficult, because of the 
combinatorial uncertainty and the fact that the top quarks from the decay of additional 
Higgs bosons with the mass of 350-500 GeV have relatively low velocity. Thus, for the mass 
determination, we may rely on the other decay modes, accompanied to the tt mode, such as 
bb and r+r“ [|68 |. 

Since the enhancement of the couplings of H and A to top quarks is true for the all four 
types of Yukawa interactions, the discrimination of the type of Yukawa interactions needs 
to see the other decay modes, such as bb for Type-H and Type-Y and for Type-H and 
Type-X. For Type-I, the other decay modes are all suppressed. Therefore, non-observation 
of the other fermionic channel may be the signature of Type-I. 

For the above two reasons, combinations of the searches of additional Higgs bosons in 
different decay modes are important. Depending on the type of Yukawa interactions and the 
value of tan/?, the decay modes into bb and r+r“ can be explored at the LHC and/or the 
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ILC m- By observing several decay modes and their event rates, tan (3 can be determined 
experimentally In the case of sin(/5 — a) < 1, the decay modes of if —)■ hh, WW, ZZ, 
and A ^ Zh can be sizable, especially for Type-I. By the combinations of the measurements 
of these modes, experimental determination of sin(/i — a) may be performed. 

Finally, we give some comments on the searches at the CLIC |]^. At the CLIC, the 
collision energy of multi-TeV, such as 3 TeV, is proposed. With i/i = 3 TeV, the mass 
reach would be extended up to about 1.5 TeV. This is totally above the scope at the LHC. 
Thus, the direct searches at the CLIC would have a great impact in any decay mode of 
additional Higgs bosons. The top quarks from the decay of such heavier Higgs bosons are 
more energetic so that the mass reconstruction by using boosted top-jet measurement can 
be realistic. In the case that H and A are produced with large velocity, the decay products 
of each Higgs boson are well separated into different hemispheres. Therefore, the mass 
of additional Higgs boson can be reconstructed by using the invariant mass of all decay 
products in one hemisphere . 


VI. CONCLUSIONS 

We have studied the direct searches of the additional Higgs bosons in multi-top-quarks 
events at the future LHC and the ILC experiments. Additional Higgs bosons are predicted 
in any kind of models with extended Higgs sectors, and the detection of them at colliders 
is a clear signature of the physics beyond the standard model. As a benchmark model of 
the extended Higgs sector, we have considered the two Higgs doublet models with discrete 
symmetry with four types of Yukawa interactions. 

At the LHC, the signals of four top-quarks events suffer from systematical uncertainties 
of estimating the SM backgrounds. Therefore, the searches can only survey the parameter 
regions with tau/d 1.5 which is disfavored by experimental constraints from flavor physics 
without adding new particle contents to the model. At the ILC, although the mass reach is 
almost limited by its beam energy, we have shown that the parameter regions with tan (3^8- 
15 can be surveyed depending on the type of Yukawa interactions. We have also discussed 
that further investigation may be performed at the ILC, such as discrimination of the type of 
Yukawa interactions and determination of the mass of additional Higgs bosons by combining 
the observations of the other decay modes of H and A in addition to the ti mode. 
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